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Abstract 



Monte Carlo (MC) simulations and series expansion (SE) data for the energy, 
specific heat, magnetization and susceptibility of the 3-state Potts model on the 
square lattice are analyzed in the vicinity of the critical point in order to estimate 
^ ■ universal combinations of critical amplitudes. We estimate these amplitudes using 

the correction-to-scaling exponents predicted by conformal field theory. We also 
form effective ratios of the observables close to the critical point and analyze how 
they approach the universal critical-amplitude ratios. In particular, using the du- 
ality relation, we show analytically that for the Potts model with q < 4 states, 
the effective ratio of the energy critical amplitudes always approaches unity linearly 
with respect to the reduced temperature. This fact leads to the prediction of rela- 
tions among the amplitudes of correction-to-scaling terms of the specific heat in the 
low- and high-temperature phases. We present numerical and analytical support 
for the form of the first two correction-to-scaling terms. Our results for the ampli- 
tude ratios closely agree with the theoretical predictions and the earlier numerical 
estimates of the specific-heat and the susceptibility amplitude-ratios. 
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1 Introduction 

The universal character of appropriate combinations of critical amplitudes is an impor- 
tant prediction of scaling theory which in some cases still remains incompletely verified 
and subject to controversies. The concept of universality is very fruitful for the classifi- 
cation of the models and of the real physical systems undergoing phase transitions. The 
set of critical exponents and critical amplitudes describes the behavior of a system in 
the vicinity of the critical point. The universal combinations of critical amplitudes, to- 
gether with the critical exponents, fully characterize the universality class [1]. The Potts 
model [2,3], one of the paradigmatic models exhibiting continuous phase transitions is 
a good frame to reconsider the question of universal combinations of amplitudes. The 
universality class of the Potts model at its critical point is determined by the number 
of states q. The two-dimensional Potts models with three and four states can be ex- 
perimentally realized as strongly chemisorbed atomic adsorbates on metallic surfaces at 
sub-monolayer concentrations [4] in accordance with the group-theoretical classification of 
the phase transitions of such systems [5-7] . Although critical exponents can be measured 
quite accurately for adsorbed sub-monolayers, confirming that these systems actually be- 
long to the three-state [8] or to the four-state Potts model classes [9], it is unlikely that 
the low temperature LEED methods can be pushed [10] to determine also the critical 
amplitudes. Therefore, the numerical analysis of these models is the only available tool 
to check analytic predictions. 

We shall restrict our analysis to the critical amplitudes and critical exponents which 
describe the behavior of the residual magnetization M, the reduced susceptibility %i an d 
the reduced specific heat C of the system at zero external field in the vicinity of the 
critical point 



M(r) w B(-rf, r < (1) 
X±(r) « r ± \r\-\ (2) 

C±(r) « — \r\~ a . (3) 
a 

Here r is the reduced temperature r = (T — T c )/T and the labels ± refer to the high- 
temperature and low-temperature sides of the critical temperature T c . In addition to 
the above quantities, for the Potts models with q > 2 a transverse susceptibility can be 
defined in the low-temperature phas^] 

Xt(t) « T t (-t)-\ (4) 

For the 2D Potts model the critical exponents have been calculated exactly [11-14] 
in terms of the number of states q. Introducing the parameter y related to the number of 
states q of the model by 

cos— = -y/q. (5) 



^dn the following we will use the notations Y l and Tt for the longitudinal and transverse susceptibility 
amplitudes in the low temperature phase. T_ also used in the literature is identified to Yl- 
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we have for the thermal scaling dimension x e = (1 — a)/v 



and for the magnetic scaling dimension x a — fl/v 

The central charge of the corresponding conformal field theory is also simply ex- 
pressed [14] in terms of y 

Analytical estimates of critical amplitude-ratios for the g-state Potts models with 
q = 1, 2, 3, and 4 were obtained by Delfino and Cardy [15]. They used the exact two- 
dimensional scattering field theory of Chim and Zamolodchikov [16] and estimated the 
ratios using a two-kink approximation for 1 < q < 3 and the contributions from both the 
two-kink approximation and the bound state for 3 < q < 4. This approximation leads 
to the value c = 0.792 for the central charge of the 3-state Potts model, for which the 
exactly known value is c = 4/5. Thus, the value of the central charge c is reproduced with 
an accuracy of one per cent. Using this approximate value in (151), we can calculate the 
scaling dimensions from ©-([7]) obtaining x a = 0.1332 and x e = 0.806, which differ from 
the exactly known dimensions (2/15 and 4/5 respectively) by less than one per cent. 

The universal susceptibility amplitude-ratios T + /Tl were calculated by Delfino and 
Cardy in [15]. Later Delfino, et al. [17] estimated analytically also the ratio of the trans- 
verse to the longitudinal susceptibility amplitude Tt/Tl. They obtained: 

r+/r L « 13.848, r T /rL~ 0.327. (9) 

In the same paper [17], MC simulations were also reported which yield 

r+/r L ^io, r T /rL = o.333(7). (10) 

In our previous paper [18], from the analysis of MC data and of old SE data for the 
3-state Potts model, we estimated T + /Tl = 14±1 in fair agreement with the theoretical 
prediction. The key point of our analysis was a fit including a correction-to-scaling term. 
Quite recently, these results were confirmed and substantially improved by Enting and 
Guttmann [19] who analyzed new longer series expansions derived by the finite lattice 
method. Their remarkably accurate estimates for the 3-state Potts model are: 

r+/r L = 13.83(9) - 13.90(15), Tt/Pl = 0.325(2) - 0.329(2) . (11) 

In the present paper, devoted to the 3-state Potts model, we present new - more 
accurate - MC data supplemented by a reanalysis of the extended series derived by Enting 
and Guttmann [19]. We also present numerical and analytical support for the form of the 
first two correction-to-scaling terms. 
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We shall be concerned with the following universal combinations of critical amplitudes 
*± I± »+-M± p- - A - V - no) 

AJ ry ry - £2 > ^ 

where the last two are a consequence of the scaling relational a = 2 — 2/? — 7. To the various 
critical amplitudes of interest, A±, T±,. . . , we have associated appropriately defined "ef- 
fective amplitudes", namely temperature-dependent quantities A±(|r|), r±(|r|), . . ., which 
take as limiting values, when \t\ — > 0, the critical amplitudes A±, T±, . . . . By analogy, 
we also considered the "effective ratios" of some amplitudes, e.g., A + (r)/A-(— r) which 
takes in the critical limit the value of the first amplitude ratio in the list (j!2p . To avoid 
any risk of confusion, reference to these temperature-dependent quantities is always made 
with their explicit r— dependence. 

A central idea in our approach is to use the duality relation in order to improve the 
estimates of the effective amplitude-ratios measured at dual temperatures. In a first anal- 
ysis, we also study directly the ratios of the effective amplitudes in (fl2j) at symmetric 
reduced temperatures ±r above and below T c . Better than separately estimating high- 
and low-temperature effective amplitudes, these effective ratios enable us to minimize 
correction-to-scaling effects in the accessible critical window. We also use the duality 
relation to estimate the leading correction-to-scaling amplitudes in the behavior of the 
specific heat and of the susceptibility. For this purpose, we compute ratios also on the du- 
ality line, e.g., the effective susceptibility amplitude ratio r +(t) /T l(t*) = X+(P) /Xl{P*) 
as the ratio of x+(P)i the high-temperature susceptibility at the inverse temperature /3, 
and Xl(P*), the low-temperature susceptibility at the dual inverse temperature (3*. 

As a final result of our analysis, we estimate the susceptibility critical amplitude-ratios 
which, for q = 3, take the values IY/IY = 0.3272(2) and T + /Tl = 13.83(8). These ratios 
are consistent with our previous results and with the predictions of Ref. [15]. 



2 Model and technical details 

The Hamiltonian of the Potts model [2] reads as 

(ij) 

where Sj is a "spin" variable taking integer values between and q — 1, and the sum is 
restricted to the nearest neighbor sites (ij) on a lattice of N sites with periodic bound- 
ary conditions. The partition function Z is defined as usual by the sum over all spin 
configurations 

conf 



2 We refer the reader to Ref. [1] for a detailed discussion of the universality of the critical amplitude 
ratios. 
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with (3 = 1/ksT, and /cb the Boltzmann constant (fixed to unity). On the square lattice 
in zero magnetic field, the model is self-dual. Denoting by f3* the dual of the inverse 
temperature @, the duality relation [2] 

(eP - 1) {e^ -l)=q (15) 

determines the critical value of the inverse temperature f3 c = ln(l + y/q). From the 
duality transformation of the partition function Z(/3) = (g _1//2 (e /3 — 1)) 2N q~ l Z '(/?*) , a 
similar transformation follows for the internal energy density E(f3) = — N -1 -^ In Z(/3). 
The values E(f3) and E((3*) of the energy density at dual temperatures are thus related 
through (see, e.g., Ref. [20]) 

(1 - e~?) E(0) + (1 - e-F) E(P*) = -2. (16) 

Dual reduced temperatures r and r* can be defined by (3 = (3 C {\ — r) and /3* = 
f3 c (l + r*). Close to the critical point, r and r* coincide through the first order, since 

T* =T+ l -^±f±T 2 + 0(T 3 ). 

2.1 Monte Carlo simulations 

We use the Wolff algorithm [21] for studying square lattices of linear size L with peri- 
odic boundary conditions. Starting from an ordered state, we let the system equilibrate 
in 10 5 — 10 6 steps measured by the number of flipped Wolff clusters. The averages are 
computed over 10 6 — 10 7 steps. The random numbers are produced by an exclusive-XOR 
combination of two shift-register generators with the taps (9689,471) and (4423,1393), 
which are known [22] to be safe for the Wolff algorithm. 

The order parameter of a microstate M(t) is evaluated during the simulations as 

qN m /N - 1 
q-l 

where N m is the number of sites % with = m at the time t of the simulation [23], and 
m E [0, 1, (q — 1)] is the spin value of the majority of the spins. iV = L 2 is the total 
number of spins. The thermal average is denoted by M = (M). 

Thus, the reduced longitudinal susceptibility in the low-temperature phase is measured 
by the fluctuation of the majority-spin orientation 

P- 1 XL = (N 2 J - (N m ) 2 (18) 

and the reduced transverse susceptibility is defined in the low-temperature phase as the 
fluctuation of the minority of the spins 

while in the high-temperature phase x+ is given by the fluctuations in all q states, 

r 1 X+ = ^E(W- W 2 )> (20) 

^ /i=0 
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where iV M is the number of sites with the spin in the state /i. Properly allowing for 
the finite-size effects, this definition of the susceptibility is, in both phases, completely 
consistent with the available series expansion data [18]. 

The internal energy density of a microstate is calculated as 



(ij) 

its ensemble average is denoted by E = (E) and the reduced specific heat per spin is given 
by the energy fluctuations, 



We have simulated the model on square lattices with linear sizes L = 20, 40, 60, 
80, 100 and 200. In each case, we have measured the physical quantities in a range of 
reduced temperatures called the "critical window" and defined as follows. Assuming a 
proportionality factor of order 1 in the definition of the correlation length, the relation 
L fa £ oc | T" | ^ yields the value of the reduced temperature at which the correlation length 
becomes comparable with the system size L and thus below which the finite-size effects 
are not negligible. This value defines the lower end of the critical window. For example in 
the q = 3 case, for systems of sizes L = 100 and L = 200, we have \r\ m i n (L = 100) ~ 0.004 
and \r\ min (L = 200) ~ 0.0017, respectively. The upper limit of the critical window is fixed 
for convenience as the value of the reduced temperature up to which the corrections to 
scaling in the Wegner asymptotic expansion [24] do not exceed a few percent, say 2 — 3%, 
of the leading critical behavior Eqs. (jT}|3]) • 

2.2 Series expansions 

Our MC study of the critical amplitudes will be supplemented by an analysis of the high- 
temperature (HT) and low-temperature (LT) expansions for q = 3 recently calculated 
through remarkably high orders by Briggs, Enting and Guttmann [19,25]. In terms of 
these series, we can compute the effective critical amplitudes for the susceptibilities, the 
specific heat and the magnetization and extrapolate them by the current resummation 
techniques, namely simple Pade approximants (PA) and differential approximants (DA) 
properly biased with the exactly known critical temperatures and critical exponents. 

The LT expansions, expressed in terms of the variable z = e~@, extend through z m in 
the case of the energy. The expansion of the longitudinal susceptibility extends through 
z 71 . In the case of the transverse susceptibility the corresponding order is z 53 . The 
magnetization expansion extends through z A7 . 

The HT expansions, computed in terms of the variable v = (1 — z)/(l + (g — l)z), 
extend to v 46 in the case of the energy, and up to v 28 in the case of the susceptibility. 

It is useful to point out that, for convenience, in Ref. [19] the product of the suscep- 
tibility by the factor q 2 /(q — 1), rather than the susceptibility itself, has been tabulated 
at HT because it has integer expansion coefficients. For the same reason, at LT the mag- 
netization times q/(q — 1) has been tabulated. Therefore the appropriate normalizations 
should be restored in order that the series yield amplitudes consistent with the MC results. 




(21) 



(3- 2 C 



BE 



((E 2 ) - (E) 2 ) N. 



(22) 
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As a general remark on our series analysis, we may point out that in the q = 3 case, 
the accuracy of the amplitude estimates is good, due to the relatively harmless nature of 
the power-like corrections to scaling. 

3 Critical amplitudes and universal ratios 

3.1 Expected temperature-dependence of the observables 

In the vicinity of the critical point, the reduced specific heat is generally expected to 
behave as 

C(r) = — |r|- Q J- corr (|r| A ) + G bt (r), (23) 
a 

where J r C orr(|'7"| A ) is the correction-to-scaling function and Quij) represents an analytic 
background (bt here stands for "background term") which accounts for non-singular con- 
tributions to the specific heat, i.e. Gu{t) = D± + -D±|t| + . . . . The specific-heat 
and the leading correction-to-scaling exponents for q = 3 are given by a = 1/3 and 
A = — v{D— x e2 ) = 2/3 where D = 2 is the space dimension [26] and x e2 = (4+2y)/(2—y) 
is the next-to- leading thermal exponent. In the HT phase and in the LT phase respec- 
tively, we can thus writ^l 

.) + D + + D' + t + ... (24) 
+ ...) +D_ +D'_\t\ + ... (25) 

where gsa,± is the amplitude of the leading correction-to-scaling, b± is the next correction 
term and so on. In the correction-to-scaling factor the ellipse denotes terms in |r| 2A , 
|r| 3A , and also terms with an increasing sequence of other exponents [14] A' = 
—u(D — x e3 ) = 5A, A" = —u(D — x €4 ) = 14A, etc. Other quantities should obey 
similar expansions including, beside the leading singularity, all corrections to scaling and 
background corrections, 

Obs.(±|r|) ~ Ampl. x \r\^ x (1 + corr. terms) + backgr. terms, (26) 
corr. terms = a|r| 2 ^ 3 + b\r\ + . . . , (27) 
backgr. terms = D + Di\r\ + . . . (28) 

where < is a critical exponent which is known and depends on the observable considered. 

3.2 Specific-heat critical amplitudes 

Although duality determines exactly the ratio of the specific-heat amplitudes A + /A- = 1, 
it is instructive to define an efficient numerical procedure to compute this ratio. Moreover, 



CAt) 



C-(~\r\) 



-^r- a (l + a A+ T A + b + r + .. 
a v 



T 



a 



(1 + a A , 



3 Here and in the following, unless the contrary is stated, r is defined as positive, but when the physical 
quantities are measured in the LT phase, r = — |r|, their temperature-dependence is explicitly denoted 
asQ(-M). 
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the actual values of the amplitudes A + and A_, albeit non-universal, are themselves 
informative, since they enter into other universal combinations, e.g. Rf = A±T±/B 2 . 

3.2.1 Corrections to scaling and background terms 

It is not convenient to extract the critical amplitudes A± directly from the specific heat. 
The energy density can be measured more accurately in the MC simulations, so in the 
following we shall study the dominant corrections to scaling and extract the background 
term from an analysis of the energy density. 

In the HT phase and in the LT phase respectively, the energy can be conveniently 
written as 



E+(t) = E + _ + r 1 "" (1 + a A , + r A + b + r +■■■)+ D + r + .. . (29) 

£_(-|t|) = E - — - . . |r| 1 - a (l + a A ,_|r| A + fe_|r| + ...) + £L|r| + . . . (30) 
a(l — a)p c v ' 

The minus sign in front of A_/a(l —a)/3 c is needed in order to recover, from the definition 
M) = j3 c dE-/dr = —f3 c dE-/d\r\, the convenient specific heat amplitude +A_/a. 
The last term represents the analytic background which may be rewritten as (A±/a(l — 
a)/? c )|r| 1_Q (i-t|r| a + . . . , in order to be incorporated in the first sum, 

E±(±\t\) = E ± A± It) 1 -" (1 + a A ,±|r| A + 6±|r| + . . . + d ± \r\ a + ...). (31) 
a(l — a)p c 

In Fig. Q] we have plotted the differences AE + = E + (t) — E and AE_ = E — E_(— |r|) 
vs. the reduced temperature r. These quantities are computed from both MC data 
(symbols) and SE (lines). By definition of the critical window, the finite-size corrections 
of the MC data are negligible in the range of reduced temperatures under study. Let us 
now define the effective amplitudes A + (t) and A_(— |r|) evaluated at symmetric reduced 
temperatures 

A+(t) = a(l-a)f3 c (E + (r)-E )r a -\ (32) 
A_(-\t\) = u{l-a)P c {E Q -E_{-\T\))\T\ a -\ (33) 

The re-summation of the series expansions of these effective amplitudes can be performed 
in various ways, all of which show a good convergence. We can compute each amplitude 
by simple PA's after performing the variable transformations w — 1 — (1 — z/z c ) A in 
the LT case (or w' = 1 — (1 — v/v c ) A in the HT case) in order to allow for the leading 
corrections to scaling. Completely consistent results are obtained by computing first-order 
inhomogeneous DA's of the amplitudes directly in the natural variables z and v. 

According to Eq. (131]) . the arithmetic mean of the effective amplitudes A+(t) and 
A-(-M) 

A{t) = \{A + {t) + A_{-\t\)) (34) 
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Figure 1: Energy differences AE + and Ai?_, calculated from MC data for system sizes L = 100 
(squares) and L = 200 (circles) and from SE data (continuous lines). 



is expected to behave as 

A(r) =A(l+ A+d+ + A A - d - r a ) + 0(t*) (35) 

where A = \{A + + A_) and the correction term comes from the leading correction to 
scaling in Eq. fl3Tj) . 

In order to compare Eq. (135]) with the numerical data, we have plotted in Figure [2] the 
effective amplitudes A + (t), A_(— |r|) vs. | t | 1//3 . The SE data for the mean A{r) are well 
represented asymptotically by the expression (solid line in Figure [2]) 

A(t) w 0.399(2) - 0.283(l)r 1/3 . (36) 

for Irj 1 / 3 > 0.16 (which corresponds to the left boundary of the critical window as dis- 
cussed above). This yields the estimate A = 0.399(2). We can also conclude from the 
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Figure 2: A + (t) = a(l - a)/3 c A J E+r a - 1 (open squares) and A_(-|r|) = a(l - a)/3 c A J E_|r| a ~ 1 
(open triangles). Their arithmetic mean A(t) as computed from SE data (open circles), or from 
MC data(closed circles), together with the fits (solid lines: fit of SE data; dotted lines: fit of 
MC data). 



essentially linear behavior of Ayr) with respect to r a that the higher-order corrections 
are rather small. Therefore, we can empirically argue that A + a& )+ ~ — A-a&-. Possibly 
also a cancellation of some higher order terms might occur in Eq. (|34"|) . This implies that 
the most important correction to be taken into account comes from the background term. 
It is necessary to include the additional terms ±0.15|r| 2 into Eq. (I5rj|) (obtained from the 
SE data) in order to catch the behavior at a larger distance from the critical point. This 
expression is more accurate also for smaller values of |r| (not accessible through MC due 
to the finite size effects). 

3.2.2 The effective amplitude-ratio A + (t)/A_(t*) on the dual line 

The previous empirical observation of a cancellation of several correction-to-scaling terms 
in appropriate combinations of effective amplitudes can be restated more rigorously by 
duality arguments. We define effective amplitudes at dual values of the reduced temper- 
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ature^, 

A+(t) = a(l-a)(3 c (E + (r)-E )T a -\ (37) 
A_(t*) = a {l- a )P c {E G -E_{r*)){r*r- 1 (38) 

and their ratio 

Mr) (E + (t)-E )t^ 

A_(t*) (E - E_{T*)){r*) a - 1 1 ' 

the constant E being the value of the energy at the transition temperature, E = E([3 C ) = 
-1 - l/y/q. 

Using the asymptotic expansions Eqs. (129]) and (!30|) in the duality equation Eq. ([TBI) 
and expanding for small r, we obtain an infinite set of relations among critical amplitudes, 
such as A + = D + = —D_, cia, + = oa,- &- = b + —2a q .. etc.. These equations 

hold for any q < 4. Therefore 



= l + (3-a)a,r + On (40) 

with a q = -E Q (5 c e-^ = =« 0.5803. 

Hereafter we shall denote by A the common value of A + and A_. 

It is interesting to check numerically the validity of Eq. ( 1401) for the asymptotic be- 
havior of the ratio of effective amplitudes evaluated at dual reduced temperatures. The 
"direct division method" of HT and LT series suggested in Ref. [27] is very effective. It 
consists in computing the quotient of the A + (v) and the A-(z) series after taking v = z 
(remember that v c = z c ). This amounts precisely to compute the ratio of the effective 
amplitudes at dual temperatures. The quotient series thus obtained is resummed by sim- 
ple PA's or DA's and can be extrapolated to the critical point obtaining the very accurate 
estimate A + /A- = 1.000000(3). The results for A + (t)/A^(t*) shown in figure [3] are com- 
pared with the same ratio computed as function of the symmetric reduced temperature 
|tU + (t)/A_(-|t|). 

The fit of MC data at dual temperatures yields a value for the slope ~ 1.5 which is 
not far from the expected value 8/3a q ~ 1.547 of Eq. (1401) . The estimate of the same 
quantity in the case of symmetric reduced temperatures \t\ yields 1.38(1). 

These plots show the cancellation of the leading non-analytic correction to scaling in 
the ratio A + /A_ of effective amplitudes. We can now take advantage of this remark to 
improve the fit of the energy data, presented in the previous subsection. First, let us 
construct the mean of the effective amplitudes evaluated at dual reduced temperatures 
Aiuaiir). For small r the SE curve has the behavior (in perfect agreement with Eq. (J36))) 

A d uai(r) = 0.399(1) - 0.283(2)r 1 / 3 (41) 



4 Notice that, with our definitions, r* is positive and characterizes the LT phase. 
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0.00 0.05 0.10 0.15 
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Figure 3: The effective-amplitude ratio A+(t)/A-(—\t\) at symmetric reduced temperatures. 
The SE data are indicated by closed circles and the MC data by open circles. The ratio 
A + {t)/A^(t*) at dual reduced temperatures is denoted by closed squares in the case of the 
SE data and by open squares for the MC data. The solid line is the prediction Eq. (|40p . 



As noticed above, the coefficient of the higher-order term at least one order of mag- 

nitude smaller. In the linear approximation of the duality relation (|40j) . A+(\t\) / A-(t*) ~ 
1 + (3 — a)a q r, one has 

Ab(±|r|) = A dua i(r)(l ± (3 - a)a q /2r). (42) 

Combining (JHJ), (j4"2"l) and using (3 — a)a q /2 0.774, one gets 

Afc(±|r|) « (0.399(1) - 0.283(2) r 1/3 )(l ± 0.774 r). (43) 

These expressions are represented in Fig. [2] by the solid lines. In order to extend the SE 
data representation of Eq. (j4ip to larger values of |r|, we can add the next background 
correction term, estimated to be =f0.15t 2 . 
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3.2.3 Energy and specific heat temperature dependence 

From Eqs. (S3), (}3"8l) . (I4T1) and ( 1421) . we get a numerical expansion of the energy of the 
3-state Potts model: 

E + (r)-E = 1 r x -M(r)(l + (3-a)a g /2r) 

a(l — ajpc 

= 1.787(5)r 2/3 (1 + 0.774(l)r - 0.412(4)r 4/3 ) - 1.269(9)r, (44) 

E -E_(-\r\) = 1 irl 1 -M(irl)(l-(3-a)a g /2lrl) 

a(i — ajpc 

= 1.787(5)|r| 2/3 (1 - 0.774(1)|^| + 0.412(4)|r| 4/3 ) - 1.269(9)|r|. (45) 

Note that the regular linear term appears through the combination of the | -7~ | 1//3 correction 
in Eq. (I4~TT) with the leading singularity in (r^ . 

As a conclusive test, we have plotted Eqs. (1441) and ( 1451) in figure [U but they cannot 
be distinguished from the solid lines representing SE data. 

We use expressions (I45ti44[) to calculate the specific heat, and plot the results in the 
figure H] together with the MC and SE data. 



3.3 Susceptibility and magnetization amplitudes 
3.3.1 The ratio T T /T L 

The transverse and longitudinal reduced susceptibilities are expected to have the following 
asymptotic form in the LT phase [24] 



|r|) 


= r T |r|-^ T (|r| A ) 


+ £r(M), 


(46) 


M) 


= r L |r|-^ L (|r| A )- 


f &(M), 


(47) 



where 7 = 13/9 [14,26] and the correction-to-scaling exponent A = 2/3 is the same as 
above. For the purpose of the fit in the low temperature range accessible by our MC 
simulation, we shall use the following expansion of the reduced susceptibility 

x (-\r\) = r|r|- 7 (l + a A \r\ A + a 2A \r\ 2A + ... + b\r\ + ... + d\r[< ...), (48) 

where for simplicity, we choose a notation for the amplitudes of the corrections to scaling 
very similar to that adopted in the previous section, since there is no risk of confusion, 
and we incorporate the leading background term whose amplitude is denoted by d inside 
the main parenthesis. The ratio of the effective amplitudes IY(— |r|) = |t| 7 Xt( — | r l) and 
= M 7 Xl(— |t|) thus behaves as 

t 

-(1 + (a A ,r - a A , L )\T\ A + (b T - b L )\r\ + 

L 

(a 2A ,T - a 2 A,L + a% L - a A , T a AiL )|r| 2A + (d T - rf L )|r| 7 + 0(|r| A+1 )). (49) 



rv(-H) 
r L (-\r\) 
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In Figure [5] we have represented the MC data for this effective ratio, plotted against 
the reduced temperature \t\. The data in Fig. [5] do not show any finite-size dependence 
even for the smallest lattice size, although their spread becomes obviously smaller for 
larger system sizes and longer MC runs. We performed a fit of the MC data collected 
for L = 100 and with the largest statistics (represented by closed circles in figure E]) to 
expression (|49|) . The values of the parameters are reported in tabled]for different trial fits. 
It is worth noticing that, while the values of the parameters are sensitive to variations of 
the limits of the critical-region window, the amplitude ratio changes only within the error 
bars. Therefore the good quality of the MC data and the large size of the critical-region 
window enable us to estimate rather accurately the amplitude and the first correction 
term in the asymptotic expansion (1491) . In table [TJ the leading correction comes from 
the regular term (last column) and the coefficient is clearly not negligible, therefore the 
fits # 2, 3 and 6 have our preference. Since it is hard to decide which one among the 
three is the most reliable, we consider that an average result Yt/Tl — 0.327(1) is a safer 
value. Figure if replotted as a function of |t| 2 / 3 shows a remarkably linear behavior for 
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Figure 5: Effective ratio of amplitudes \t\)/Tl(— |t|) for lattices of linear sizes L = 20 

(boxes), L = 40 (up triangles), L = 60 (down triangles), L = 80 (diamonds), L = 100 (stars) 
computed with Nmc = 10 5 Monte Carlo steps. We have indicated by closed circles the results 
of a computation on a L = 100 lattice with Nmc = 10 6 Monte Carlo steps. The dashed line 
(hardly visible) represents a fit of the data for L = 100 to Eq. (|49p . The values of the coefficients 
are listed in the last entry of table [TJ The solid line represents the [22,22] Pade approximant to 
the ratio of the LT series in a conveniently transformed variable. 



|r| 2/3 < 0.2, while the result of the linear fit: T T /T L = 0.329(1) supports our choice of 
the value 2/3 for the leading confluent correction exponent. 

Although possibly less accurate, the MC data for the low temperature susceptibilities 
can be fitted in the same window, leading to 

Xt = 0.00409(2)|t|~ 13/9 (1 + 0.32(7)|r| 2/3 - 1.88(ll)|r|), (50) 
Xl = 0.01270(4)|r|- 13/9 (l - 1.60(5)|r| 2/3 + 0.44(8)|r|). (51) 

The ratio of the amplitudes, T T /T L w 0.322, is fairly consistent with the value quoted 
above, but since no background term is included we consider that this value is not highly 
reliable. The differences between coefficients a^'s and 6's appearing in Eq. (!50j) and (!5T!) 
are also compatible with those quoted in table [TJ 

The quotient of the LT series for Ft and Tl can also be studied either most simply by 
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Table 1: Results of the fit to the MC data for the ratio Tt/Tl of the transverse and the 
longitudinal susceptibility in the critical region window. The absence of an entry in the table 
indicates that we have not included the corresponding parameter in the fit procedure. 



Fit # 


ratio 




corrections coefficients 


r T /r L 


<*M A 


tx M 


oc |rpA 


oc |t|t 


1 


0.328(2) 


1.04(26) 


1.11(87) 


-2.03(79) 




2 


0.328(2) 


1.13(23) 


0.56(6) 




-1.63(65) 


3 


0.327(2) 


1.24(12) 




1.9 ±2.6 


-3.1 ± 2.7 


i 


0.324(1) 


1.71(5) 


-1.10(7) 






5 


0.326(1) 


1.37(2) 




-1.02(6) 




6 


0.326(1) 


1.33(2) 






-1.08(6) 



PA's in the transformed variable w as mentioned above or by DA's in the variable z. Both 
kinds of approximants are smoothly extrapolated to the value IY/T^ = 0.3272(2) for the 
critical amplitude ratio. The critical amplitudes can also be separately computed by DA's 
obtaining the estimates IV = 0.004166(5) and = 0.01273(1). Our results, summarized 
and compared with those of previous studies in Table [2] and [31 agree completely with the 
prediction T T /T L = 0.327 by Delfino, et al. [17] whose uncertainty is presumably of the 
order of a half percent. 

3.3.2 The critical amplitude of the magnetization. 

The magnetization is expected to behave in the critical region as 

M(-\t\) = B\rf (1 + a A \r\ A + ... + b\r\ + ...)+ D\r\ + ..., (52) 

where /3 = 1/9 is the magnetization exponent and as usually we have indicated only the 
leading corrections to scaling and the analytic background. 

The analysis of the LT expansion for the effective amplitude B(— |r|) = |r|~^M(— |r|) 
of the magnetization can be efficiently performed either by simple PA's in the transformed 
variable w or by DA's. We obtain the estimate B = 0.819(1) for the critical magnetization 
amplitude. DA's also indicate that the exponent of the leading correction to scaling is 
0.73(3), 

B(-\t\) = 0.819(1) - 0.226(2)|r| a73(3) . (53) 

Comparing (152"]) to (1531) one can identify | -7- 1 o.T3(3) w jth |r| A . As an illustration, a fit (with 
the exponent fixed to the value A = 2/3) of the SE data is shown in the figure [6] (dotted 
line) and compared with the fit fl53|) (solid line). 

A highly compatible fit is obtained for the MC data which yields 

M = 0.818(l)|r| 1/9 (l - 0.65(3)|r| 2/3 - 0.400(7)|r| + 0.2865|r| 4/3 ). (54) 
3.4 ratios 

We can finally estimate the value of the universal amplitude ratio Rq = A + T + /B 2 . If 
we use our estimates A + = 0.399(1), B = 0.819(1) and the value T + = 0.176(1) obtained 
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\T\ 

Figure 6: The effective amplitude of the magnetization M as computed from SE (open circles) 
together with the fits (lines, see in text). Insert: The magnetization M as function of It) 1 / 9 
(open circles) and a fit to the curves ^^^(l + a>M\T\ Actt ) with the exponent A e g fixed to the 
values A e ff = A = 2/3 (upper line) and A e g = 1- /3 = 8/9 (lower line). 



from the presently available series [19] (much longer than those used in our previous 
study [18]), we conclude 

R£ = 0.1049(8). 

This value compares very well with the estimate = 0.1041 by Delfino and Cardy [17]. 

The corresponding Monte Carlo estimates are A + = 0.396(9), B = 0.818(1) and 
T + = 0.1783(7) leading to a value for the ratio i?+ = 0.1054(29) consistent with SE 
estimation albeit less accurate. 

An alternative approach, leading to a very similar numerical estimate R^ = 0.1043(8), 
consists in expressing Rc as the combination of the effective amplitudes A(v), T(v) and 
B(z), after taking v — z, namely on the dual line in the z — v plane. This procedure, 
however, wastes a large part of the available series coefficients. 

Using the value T L = 0.01273(1) obtained from SE one can estimate R£ = 0.00756(4) 
while using the MC value T L = 0.01270(4) one obtains R^ = 0.00751(18). These values 
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agree well with Delfino and Cardy [17] analytic prediction R c = 0.00752. 

3.5 Summary of the results for the 3-state Potts model 

By MC simulations and by series extrapolations, we have computed all universal amplitude- 
ratios studied in Refs. [15] and [17] for the q = 3 Potts model, with the exception of those 
involving the correlation length. We have shown that the cancellation of the leading 
non-analytic corrections to scaling in the ratio of the effective amplitudes of the energy 
evaluated at dual temperatures leads to a very accurate estimate of this amplitude ratio. 
Using this result, we have given an asymptotic numerical representation of the energy in 
a vicinity of the critical temperature, and observed that the main correction to scaling is 
due to the background term. 

In table [2] we have collected the results from our fits of the main observables to the 
generic form given in Eq. (1261) . These fits are performed either using MC data or SE 
data, and their results are reported in the last column. From the estimates obtained, 
we can form the universal combinations reported in table [31 We have in general found 
more accurate results by fitting effective amplitude ratios rather than the amplitudes 
themselves. A comparison of our numerical results with the analytical predictions and 
with numerical estimates from other sources generally shows a very good agreement. 

As a conclusion, we would like to mention that our primary interest is also to study 
the universal combinations of amplitudes in the case of the 4— state Potts model for 
which the available results are less convincing due to the presence of logarithmic correc- 
tions to scaling. We believe that the analysis proposed here in the case of the 3— state 
model (essentially based on an analysis of quantities at dual temperatures and of effective 
temperature-dependent combinations) can be successfully transposed to the 4— state case 
(see Refs. [28,29] for preliminary results). 
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